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1. Introduction. During the last decades, a wealth of analytical and numerical 
tools from control theory and optimization have been successfully employed to ana- 
lyze and control the performance of quantum mechanical systems, advancing quantum 
technology in areas as diverse as physical chemistry, metrology, and quantum infor- 
mation processing [T]. Although measurement- based feedback control [5] and the 
promising coherent feedback control [3] have gained considerable attention, open-loop 
control has been proven quite effective. Controllability results for finite- and infinite- 
dimensional quantum mechanical systems have been obtained, clarifying the control 
limits on these systems [iHllj. Analytical solutions for optimal control problems de- 
fined on low-dimensional quantum systems have been derived, leading to novel pulse 
sequences with unexpected gains compared with those traditionally used [T^Hin] • And 
numerical optimization methods, based on gradient algorithms or direct approaches, 
have been used to address more complex tasks and to minimize the effect of the 
ubiquitous experimental imperfections |21H29j . 

At the heart of modern quantum technology lies the efficient cooling of trapped 
atoms, since it has created the ultimate physical systems thus far for precision spec- 
troscopy, frequency standards, and even tests of fundamental physics [30], as well as 
candidate systems for quantum information processing |31) . In the present article 
we study a time-optimal control problem related to the frictionless cooling of atoms 
trapped in a time-dependent harmonic potential. Frictionless atom cooling in a har- 
monic trapping potential is defined as the problem of changing the harmonic frequency 
of the trap to some lower final value, while keeping the populations of the initial and 
final levels invariant, thus without generating friction and heating. Conventionally, 
an adiabatic process is used where the frequency is changed slowly and the system 
follows the instantaneous eigenvalues and eigenstates of the time-dependent Hamilto- 
nian. The drawback of this method is the long necessary times which may render it 
impractical. A way to bypass this problem is to use the theory of the time-dependent 
quantum harmonic oscillator |32j to prepare the same final states and energies as with 
the adiabatic process at a given final time, without necessarily following the instan- 
taneous eigenstates at each moment. Achieving this goal in minimum time has many 
important potential applications. For example, it can be used to reach extremely low 
temperatures inaccessible by standard cooling techniques [33], to reduce the velocity 
dispersion and collisional shifts for spectroscopy and atomic clocks [34J, and in adia- 
batic quantum computation i35l . It is also closely related to the problem of moving in 
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minimum time a system between two thermal states, as for example in the transition 
from graphite to diamond |36j . 

It was initially proved that minimum transfer time for the aforementioned problem 
can be achieved with "bang-bang" real frequency controls [36 . Later, it was shown 
that when the restriction for real frequencies is relaxed, allowing the trap to become 
an expulsive parabolic potential at some time intervals, shorter transfer times can 
be obtained, leading to a "shortcut to adiabaticity" [37]. In our recent work [38] . 
we formulated frictionless atom cooling as a minimum-time optimal control problem, 
permitting the frequency to take real and imaginary values in specified ranges. We 
showed that the optimal solution has again a "bang-bang" form and used this fact to 
obtain estimates of the minimum transfer times for various numbers of switchings. In 
the present article we complete our previous work by fully solving the corresponding 
time-optimal control problem and obtaining the optimal synthesis. As the terminal 
point in the problem is varied, a rather unconventional and interesting switching 
structure involving cut-loci and discontinuous switching curves is revealed. 

2. Formulation of the problem in terms of optimal control. The evolu- 
tion of the wavefunction ipit^x) of a particle in a one-dimensional parabolic trapping 
potential with time- varying frequency u){t) is given by the Schrodinger equation 



dtp 



2^dx^ ^ 2 ^ 



tP, (2.1) 



where m is the particle mass and h is Planck's constant; a; G R and tp is sl square- 
integrable function on the real line. When uj{t) is constant, the above equation can 
be solved by separation of variables and the solution is 

oo 

V'(t,x) = ^c„e-^<*/''vI;-(^), (2.2) 



where 



are the eigenvalues and 



K= + n = 0,l,... (2.3) 



are the eigenfunctions of the corresponding time-independent equation 



2 



2m dx^ ^ 2 ' " " 



Here Hn in (|2.4p is the Hermite polynomial of degree n. The coefficients c„ in 
can be found from the initial condition 

tP{0,x)'i/';;^{x)dx. 

Consider now the case shown in Fig. 12.11 where Ld{t) — for i < and uj{t) = 
ujt < for t > T. This corresponds to a temperature reduction by a factor w^/wo. 



TIME-OPTIMAL FRICTIONLESS ATOM COOLING 



3 




T 



Fig. 2.1. Time evolution of the harmonic trap frequency. 

if the initial and final states are canonical [37]. For frictionless cooling, the path uj{t) 
between these two values should be chosen so that the populations of all the oscillator 
levels n = 0,l,2,...fort>T are equal to the ones at t = 0. In other words, if 

oo 

V'(0,a;) = ^c„(0)*^°(^), 

and 

oo 

V'(t,a;) = ^c„(i)*r(a;),t>T, 

n=0 

then frictionless cooling is achieved when 

|c„(t)|2 = |c„(0)p, i >T,n = 0,l,2,... (2.5) 

Among all the paths uj{t) that result in (I2.5p . we would like to find one that achieves 
frictionless cooling in minimum time T. In the following we provide a sufficient con- 
dition on uj{t) for frictionless cooling and we use it to formulate the corresponding 
time-optimal control problem. 

Proposition 2.1. If uj{t), with uj{0) = wo and uj{t) = w(T) ^ ojt for t > T is 
such that the Ermakov equation 139}/ 

bit)+u:'mt)--^^ (2.6) 

has a solution b{t) with 6(0) = 1,6(0) = and b(t) = b(T) = (loq/ujt)'^/'^ ,t > T, then 
condition \2. 5\] for frictionless cooling is satisfied. 

Proof. Without loss of generality we assume that the initial state is the eigcnfunc- 
tion corresponding to the n-th level -0(0, x) = ^^"(x). We will show that when the 
hypotheses of PropositionlOhold then ^i){t,x) = e*""(*)5'5^^(a;), t > T, where a„(t) is 
a global (independent of the spatial coordinate x) phase factor. This and the linearity 
of (EH) imply that if 0(0, x) = J^n^o Cn(0)*;;^« (x) then V(i, x) = J^Zo c„(0)e*""(*) x 
'^n'^{x),t > T, thus condition is satisfied. 
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The frequency variations in the trapping potential change the time and distance 
scales and motivate the use of the following "ansatz" , introduced by Kagan et al. [30] , 
in (EH 



(r, x) exp 



.mx"^ b{t) 



where x = x/b{t), r = T{t) is a time rescaling, and the distance scale b{t) satisfies 
(I2.6|) and the accompanying boundary conditions. We obtain 



2™ dx^ 

If we choose the time scale T{t) such that 

dt' 



m{b + uj'^b)b^ 2 



T{t) 



b^it'Y 



(2.7) 



(2.8) 



then (|2.7p becomes 



OT 



2m c^x^ 2 

with the initial condition 0(0, x) = (x)- So (/"(r, x) = e~'-^"°^/'"'*^» (x) and 



V'(i, a;) = exp 



.mx^ b{t) 



X exp 



We will show that for t>T, where b{t) = (wo/wt)^^^, fpitjx) has the desired form. 
We examine separately each of the three terms in (|2.9p . Since b{t) = in this time 
interval, the first exponential is equal to unity. About the second exponential, observe 
from dSIH]) that 

rit)^riT) + '^it-T), 
since b{t) = {loq/cut)'^^'^ ,t > T. Also, from we have E';l° = E^-^ujo/ojt- Thus 



-iE^OT{t)/h 



g-j_E7 T{T)/h^-iE^'^ {t-T)/h 



The last term in 



satisfies 



UJQ 



1/4 



OJQ 



as it can be verified using (|2.4p . Putting all these together we see that il;{t, x) has the 
desired form for t > T. □ 

In order to find the path io{t),0 < t < T, that accomplishes frictionless cooling in 
minimum time T, we express the problem using the language of optimal control, incor- 
porating possible restrictions on u){t) due, for example, to experimental limitations. 
If we set 



b , , cj2(t) 
xi = 6, a;2 = — , u[t) = ^, 



(2.10) 
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and rescale time according to tnew = '^o^old' obtain the following system of first 
order differential equations, equivalent to the Ermakov equation (j2.6p 

ii=X2, (2.11) 

X2 = —uxi H — 3 . (2-12) 
^1 

If we set 7 = (wq/wt)^^^ > 1, the time optimal problem takes the following form 

PROBLEM 1. Find —ui < u{t) < U2 with u(0) = \,u{T) = 1/7^ such that starting 
from (a;i(0), 0:2(0)) = (1,0), the above system reaches the final point {xi{T),X2{T)) = 
(7,0), 7 > 1, in minimum time T. 

The boundary conditions on the state variables (xi, X2) are equivalent to those for 
6, 6, while the boundary conditions on the control variable u are equivalent to those for 
uj, so the requirements of Proposition 12.11 are satisfied. Parameters ui,U2 > define 
the allowable values of u{t) and it is U2 > u{0) = 1. Note that the possibility w^(t) < 
(expulsive parabolic potential) for some time intervals is permitted, Chen et al. |37) . 
It is natural to consider that also ui > 1, i.e. we can at least achieve the negative 
potential V{x) = —mujQx'^/2. Finally observe that the above system describes the 
one-dimensional Newtonian motion of a unit-mass particle, with position coordinate 
xi and velocity X2. The acceleration (force) acting on the particle is —uxi + l/a^i- 
This point of view can provide useful intuition about the time-optimal solution, as we 
will see later. 

In the next section we solve the following optimal control problem 
PROBLEM 2. Find —ui < u(t) < U2, with ui,U2 > 1, such that starting from 
(a;i(0), X2(0)) = (1,0), the system above reaches the final point {xi(T),X2(T)) = 
(7,0), 7 > 1, in minimum time T. 

In both problems the class of admissible controls formally are Lebesgue measur- 
able functions that take values in the control set [— ui,M2] almost everywhere. How- 
ever, as we shall see, optimal controls are piecewise continuous, in fact bang-bang. 
The optimal control found for problem [2] is also optimal for problem [1] with the ad- 
dition of instantaneous jumps at the initial and final points, so that the boundary 
conditions u{0) = 1 and u{T) = 1/7'* are satisfied. Note that in connection with Fig. 
12. 1[ a natural way to think about these conditions is that u{t) = 1 for t < and 
u{t) = 1/7'* for i > T; in the interval (0,T) we pick the control that achieves the 
desired transfer in minimum time. 

3. Optimal Solution. The system described by (|2.1ip . (I2.12p can be expressed 
in compact form as 

i = f{x) +ug{x), (3.1) 
where the vector fields are given by 









( 1/^? ) 




. -2^1 ) 



(3.2) 



and X E V = {{xi,X2) € : a;i > 0} and u E U = [—7/1,^2]. Admissible controls 
are Lebesgue measurable functions that take values in the control set U. Given an 
admissible control u defined over an interval [0, T], the solution x of the system p.ip 
corresponding to the control u is called the corresponding trajectory and we call 
the pair (a;, u) a controlled trajectory. Note that the domain V is invariant in the 
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sense that trajectories cannot leave T). Starting with any positive initial condition 
a;i(0) > 0, and using any admissible control u, as xi 0+, the "repulsive force" l/xf 
leads to an increase in xi that will keep xi positive (as long as the solutions exist). 

For a constant Aq and a row vector A = (Ai,A2) G (M^) define the control 
Haniiltonian as 

H = H{Xo,X,x,u) = Ao + (A, /(x) +ug{x)). 

Then the conditions of the Pontryagin Maximum Principle Jll provide the following 
necessary conditions for optimality: 

Theorem 3.1 (Maximum principle for control affine time-optimal problems). 
141^ Let {xt:(t),u^{t)) be a time-optimal controlled trajectory that transfers the initial 
condition x{0) — Xq into the terminal state x{T) = xt- Then it is a necessary 
condition for optimality that there exists a constant Aq < and nonzero, absolutely 
continuous row vector function X(t) such that: 

1. A satisfies the so-called adjoint equation 

OH 

Kt) = ~-^{XQ,Xit),x4t),u,{t)) = -{X{t),Df{x4t))+u,it)Dg{x4t))) 

2. For < t < T the function u H- _ff (Aq, A(t), x*(t), u) attains its maximum 
over the control set U at u = Ut:{t). 

3. H{Xo,X{t),x4t),u,it)) = 0. 

We call a controlled trajectory {x, u) for which there exist multipliers Aq and X{t) 
such that these conditions are satisfied an extremal. Extremals for which Aq = are 
called abnormal. If Ao < 0, then without loss of generality we may rescale the A's and 
set Ao = —1. Such an extremal is called normal. Abnormal extremals typically corre- 
spond to some degeneracies in the structure of the optimal solution (often the value 
function is no longer differentiable along these paths), but they cannot be excluded 
a priori for time-optimal control problems. For example, the solution to the time- 
optimal control problem to the origin for the harmonic oscillator, a simple text book 
example, is largely characterised by two optimal abnormal controlled trajectories. 

For the system (|2?TT1) . (f2T2|) we have 

iJ (Ao, A, X, u) — Xo + X1X2 + X2 ( — — xiuj , (3.3) 



1 



and thus 

A = -A 



1 A , / 



^ +"l -1 



Observe that H is a linear function of the bounded control variable u. The 
coefficient at u in _ff is — A2a::i and, since xi > 0, its sign is determined by $ = — A2, 
the so-called switching function. According to the maximum principle, point 2 above, 
the optimal control is given hy u — —ui if $ < and by u = M2 if $ > 0. The 
maximum principle provides a priori no information about the control at times t 
when the switching function $ vanishes. However, if ^{t) = and <b{t) ^ 0, then 
at time t the control switches between its boundary values and we call this a bang- 
bang switch. If $ were to vanish identically over some open time interval / the 
corresponding control is called singular. 

Proposition 3.2. For Problem\^ optimal controls are bang-bang. 
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Proof. Whenever the switching function $(t) — —X2{t) vanishes at some time 
t, then it follows from the non-triviality of the multiplier X{t) that its derivative 
<^{t) = —X2{t) = Xi{t) is non-zero. Hence the switching function changes sign and 
there is a bang-bang switch at time t. □ 

Thus optimal controls alternate between the boundary values u = ~ui and u — U2 
of the control set and we shall see below that the number of switchings remains 
bounded on compact subsets of the domain T). Chattering controls that would have 
infinitely many switchings on a finite interval are not possible. 

Definition 3.3. We denote the vector fields corresponding to the constant hang 
controls ~ui and U2 by X = f — uig and Y ~ f -\- U2g, respectively, and call the 
trajectories corresponding to the constant controls u = —ui and u = U2 X- and Y- 
trajectories. A concatenation of an X -trajectory followed by a Y -trajectory is denoted 
by XY while the concatenation in the inverse order is denoted by YX . 

In this paper we establish the precise concatenation sequences for optimal controls 
and in particular calculate the times between switchings explicitly. 

Proposition 3.4. All the extremals are normal. 

Proof. If (x, u) is an abnormal extremal trajectory that has a switching at time 
t, then, since X2{t) — 0, it follows from H — that we must have X2{t) — 0. The 
starting point is (1,0) and suppose that u — — ui initially. From (|2.12p it is ±2 > 
so a;2 > and a switching at a point with X2(t) > 0, not allowed for an abnormal 
extremal, is necessary in order to reach the target point (7,0). liu = U2 initially, then 
2^2(0) = 1 — U2 < and X2 < for some time interval. During this time it is ii < and 
consequently xi < 1 < 7. A switching is necessary, which takes place on the xi-axis 
for an abnormal extremal. The control changes to u = —ui and the situation is as 
before, where one more switching is necessary at a point with X2{t) > 0, forbidden for 
abnormal extremals. Thus, there are no abnormal extremals in the optimal solutions. 
□ 

We henceforth only consider normal trajectories and set Aq = —1. For normal 
extremals, H — then implies that for any switching time t we must have Xi{t)x2{t) = 
1. For an XY junction we have <i>(i) = Xi{t) > and thus necessarily X2{t) > and 
analogously optimal YX junctions need to lie in {x2 < 0}. We now develop the 
precise structure of the switchings in a series of Lemmas. We start with computing 
the evolution of the state xi{t) along an X- or F-trajectory. 

Lemma 3.5 (Time evolution of xi). The time evolution of Xi along an X- 
trajectory in the upper quadrant starting from (a, 0) is 

Mt) = Jl (a^ - +Uc^' + cosh(2ViZri), (3.5) 

while the corresponding evolution along a Y -trajectory in the lower quadrant starting 
from {l3, 0) is 



x.it) . (;32 + _L_) + 1 (;32 _ _L_) eos(2V^t), (3.6) 
Proof. A first integral of the motion along the X-trajectory is 



2 _ 2 I _ 

Xn U\X-t \ Q C, 

XI 



(3.7) 
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where c — ^uic? -^Xjo? . From (j2.12p we observe that ±2 is positive for u = —U\ and 
since a;2(0) = it follows that X2(t) itself is positive. Hence 

\/uixf + cx? — 1 
X2 = 

Xl 

and (|2.1ip gives 



\/uix'i + cx'i — 1 
Xl = — . 

Xl 

Making a change of variables according to 



the previous equation becomes 

dy 



2uixl+c 

, , (3-8) 



Integrating and using 7/(0) — 1 we obtain that 



2Juldt. 



ln(y + vV^) = 2^it. 



From this and p.Sp . equation p.Sp easily follows. 

Similarly, a first integral of the motion along the ^-trajectory is given by 

x\ + U2x\ H — n = c, (3.9) 

XI 

where now c — U20^ + We are interested in the part of the trajectory in the 

lower quadrant, X2 < 0, and thus 

\/-U2xl + cxl -1 

X2= 

Xl 

and 



\/-U2x\ + CX\-1 

Xl — . 

Xl 

If we now make the change of variables 



Vc^ - 4U2 ' 
we obtain 

— , ^ — —2\/uidt. 

Integrating this and using y{0) = 1 wc find that 

y = cos{2y/u^t). 



2u2x\-c 

y = (3.10) 
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From this and p.lOp we can easily derive (|3.6I) . Note that in the calculation of y{0) 
we used that for evolution in the lower quadrant it necessarily holds that i;2(0) < 
U2l3^ > 1/13^ so Vc^ - 4u2 = "2/3^ - l//3^ □ 

The times between consecutive switchings along optimal controls are determined 
by specific relations that we now derive. 

Lemma 3.6 (Inter-switching time). Let p = {xi,X2) be a switching point and t 
denote the time to reach the next switching point q. Ifpqisa Y -trajectory, then 

cos(2ViIJr) = ^^^^ (3.11) 



■ u^o / — ^ 2^/0^x1x2 , — , xl+uixl 
smh(2yVrr) = g j, cosh(2-y/uiT) = — ^. (3.12) 

Note that the inter-switching times depend only on the ratio X2/X1. 

Proof. These formulas are obtained as an application of the concept of a "conju- 
gate point" for bang-bang controls as originally defined by Sussmann in [32] and [13] . 
For additional background on the synthesis of optimal controlled trajectories in the 
plane, we also refer the reader to the monograph 04] by Boscain and Piccoli that 
gives a comprehensive introduction to the theory of optimal control for 2-dimensional 
systems. In an effort to make the paper self-contained, we include Sussmann's argu- 
ment. 

Without loss of generality assume that the trajectory passes through p at time 
and is at q at time t. Since p and q are switching points, the corresponding 
multipliers vanish against the control vector field g at those points, i.e., (A(0), g{p)) = 
{X(t), g{q)) = 0. We need to compute what the relation {\(t), g{q)) — implies at 
time 0. In order to do so, we move the vector g{q) along the K-trajectory backward 
from q to p. This is done by means of the solution w(t) of the variational equation 
along the y-trajectory with terminal condition w{t) = g{q) at time r. Recall that 
the variational equation along Y is the linear system w — Aw where A is given in 
(13. 4|) . Symbolically, if we denote by e*^(p) the value of the F-trajectory at time t 
that starts at the point p at time and by (e~*^)* the backward evolution under the 
linear differential equation w — Aw, then we can represent this solution in the form 

w{Q) = {e-^'')Mr) = {e-^^)*9{q) = {e^^^Uie^'^ {p)) = (e"^^). o g o e-^(p). 

Since the "adjoint equation" of the Maximum Principle is precisely the adjoint equa- 
tion to the variational equation, it follows that the function t M- {X{t),w{t)) is constant 
along the F-trajectory. Hence {X{t), g{q)) = implies that 

(A(0), ^0)) = (A(0), (e-^ag(e-^(p))) = 

as well. But the non-zero multiplier A(0) can only be orthogonal to both g{p) and 
w{0) if these vectors are parallel, g{p)\\w{0) — (e~'^^)*(jf(e'^^(p)). It is this relation 
that defines the switching time. 

It remains to compute w{0). For this we make use of the well-known relation [45] 

(e-"^), o.goe"^ = e""'^^(g) (3.13) 

where the operator adY is defined as adY{g) — [Y, g], with [, ] denoting the Lie bracket 
of the vector fields Y and g. This representation is a consequence of the fact that 



■ fr, / — ^ 2^x1x2 

x^ + U2XI 

while, ifp^ is an X -trajectory, then 
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the derivative of the function x '■ t ^ (e at t = is given by \Y,g\{p) 

and iteratively the higher order derivatives of x at are given by x''"''(0) = o,d"-Y{g) 
where, inductively, a(F''Y{g) = \Y,a(F-~^Y{g)\. For our system, the Lie algebra C 
generated by the fields / and g actually is finite dimensional: we have 

= („^; 

and the relations 

[/,[/, .9]] = 2/, [<?,[/, 5]] --25 

can be directly verified. Using these relations and the analyticity of the system, 
^tadY can be calculated in closed form from the expansion 

e"'''''{g) = J2-ad^Y{g). (3.14) 

It is not hard to show that for n = 0, 1, 2, . . ., we have that 

ad2"+ir(.g)-(-4«2)"[/,.9] 

and 

ad^-+'Y{g)^2{-4u2nf-u,g), 

SO that 

By summing the series appropriately we obtain 



5* "'''(5) =9+ sin(2V^i)[/,.9] + -^[1 - cos(2V^<)](/ - U2g). 

20i2 2m2 



Hence the field w{0) = (e ^^)*g{e'^'^ (p)) is parallel to g{p) ~ (0, —xi)'^ if and only if 
^/u^xi sin(2Y^r) + X2[l — cos(2Y^r)] = 0. 

Hence 

sin(20i;T) = ^^[1 - cos(2V^r)] (3.15) 



from which p. lip follows. Note that the solution cos(2y^T) = 1 is rejected because 
it corresponds to r = or r = "^/y^, the latter being the period of the closed 
trajectory. 

In the case of an X-trajectory the corresponding inductive relations are 

ad2"+ix(g) = (4wi)"[/,ff] 

and 

ad'"+^X{g)=2i4u,nf + uig) 
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Fig. 3.1. Consecutive switching points lie on two opposite-slope lines through the origin. Blue 
curves correspond to X -segments, red curves to Y -segments. 



for n = 0,1,2,. 



and 



A adX 



(g) =9 + sinh(2V^i)[/,5] + -L[cosh(2Viirt) -!](/ + wi.g). 

Z^ui 2ui 



For t — T this field is parallel to g at p if and only if 

^/uiXi sinli(2y^T) + a;2[cosli(2y'?IiT) — 1] = 0, 
from which we find 



sinh(2Y^T) = — 



X2 



UiXi 



-[cosh(27^T) - 1]. 



(3.16) 



Using this relation we obtain (|3.12p . The solution cosh(2Y^T) = 1 corresponds to 
T = and is rejected. □ 

Lemma 3.7 (Main technical point). The ratio of the coordinates of consecutive 
switching points has constant magnitude but alternating sign, while these points are 
not symmetric with respect to the xi-axis. 

Proof. Consider the trajectory shown in Fig. l3.1l with switching points (k, //), (C, ^) 
and (A, v). Note that (k, ji) is the intersection of an X-trajectory, passing from (a, 0) 
and a F-trajectory passing from (/3,0). From the first integrals p.7p and (13.91) we 
have 



2 2 ^ 
— UlK H ^ 



1 



-uia 



+ U2k'^ 



From these equations we obtain 



/i ^(^2 — a^){uia^K'^ + 1) 



and 



(ui + U2)a K — uia + 1 



(3.17) 
(3.18) 

(3.19) 
(3.20) 
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Note that /i > since this kind of switching can occur only in the upper quadrant. 
We win show that — —fi/n. Starting from {K,fJ,), let tq,Ts denote the time to 
reach the points (/3, 0), {(, respectively. Observe that C satisfies (|3.6|) for t — Ts~tq 
while K satisfies this equation for t = — tq- From the first relation we get 



2 V" ' 

I (f^^ ^ ) [cos(2V^ro) cos(2^T,) + sin(2V^To) sin(2V^T,)], (3.21) 

while from the corresponding relation for k 



'^^ 4 (^^ + ^) + H^' ~ ^ ^ cos(2V^.o). 



Using p.20p we obtain 



1 / „9 1 \ , , , (lio — Ml )Q!^K^ + Ml — 1 



and after an elementary but a bit lengthy calculation 



, f - sm(2V^.o) = L >. (3.23) 



Here we used that U2P^ > l//3^, since ±2 < at (/3,0), and sin(2Y^ro) > 0, since 
To < Tq/2, Tq = I ^Ju^ being the period of the closed orbit. For the terms involving 
the switching time we use p. lip . p.l7p to obtain 



COS 2yiIjT, = ( ^ \ 2 ^^^^ 2 3.24 



and 



2ie-^U2a^{n'^ - d^){uxa^K^ + 1) 
sin(2Vw2Ts) = -- — _J N 2 4 1/1 4w 2 ■ (3.25) 



Observe from (jXTTI) that sin(2^rs) < since /x > 0. By using ([3221), dSSS]), (0211), 
(P:^ in dnSH), and the relations ((X^ and 



?'+W2C' + «^^2/3' + ^, (3.26) 

C = , (3.27) 

y (ui + u-2)cP-i^'^ + (1 ~ uia*)^^ — 

^ ^ V(^^ -a2)(^^a2^2_^X) 

(wi + M2)a^K^ + (1 — U\OL^)k?- — O? 



we obtain 



so 



^ -^Z (k^ — 0?^(\l\0? + 1) 



(3.29) 
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Obviously, it is C 7^ « in general so (C, 7^ ('^j "A*): the subsequent switching point 
is different from the symmetric image of the previous switching point with respect to 
the xi-axis. 

A similar computation shows that v/X — —^/C- Note that (C,C) belongs to the 
X-trajectory passing from (ai,0), so 

e-uiC'' + i=-uial + \. (3.30) 

Using (|3:26)) . t^M\ we obtain 

2 1 (m + U2)f3H^ ~ U2/3^ - 1 
"1 2= m. • ^-^^ 

and an alternative expression for ^/C 



Starting from (C,C) l^t to,Ts now denote the time to reach the points (q:i,0), (A, z^), 
respectively. Observe that A satisfies p.Sp for t = — tq while satisfies this equation 
for t = —tq. From the first relation we get 

\2 1 ^ 2 1 
= 7: "l 9 

2 V uia? 

1/^2. 1 
- I a. ' 



2 V"' Mia? 



[cosh(2y^To) cosh(2Y^rs) — sinh(20Iiro) sinh(2Y^Ts)], 

(3.33) 



while from the corresponding relation for and by using p.3ip we obtain 

1^2, 1 \ u/o ^ N (mi - U2)^'C' + M2/3^ + 1 



2 



cosh(2V^ITTo) = n\2 ' (3-34) 



^ J -h(2V^ro) = W -;\^2- ^- (3.35) 



For the terms involving the switching time we use p.l2p . p.26p and find 

uo r- A (mi-m2)/32C" + (1 + m2/3^)C^-/32 ,„ 

cosh 2/^r, = — — - — .r,2M^n ^ ^4^72 ^' 3.36 

-(Wl + M2)P C + (1 + M2P )C - P 



smh(2ViirT,) = ... , ^.2.4^n ^„.«4V2_/^2 - (3-37) 



-(ui+M2)/32C^ + (l+M2/34)C2-/3 

By using ([g^i]) . (P35|) . ([535)1 . (jXTfl) in (1535)) . and the relations (|53T)) and 



-uiX^ + ^ = -uial + \, (3.38) 
A^ af 



we obtain 



/3C 



^-(Ul + U2)(3^C^ + (1 + M2/?4)C2 - /32 

^iP^~e)iu2P^e~i) 

Cv/-(mi + M2)/32C4 + (1 + U2P^)e - 



(3.39) 
(3.40) 
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Fig. 3.2. Blue curves correspond to X-segments, red curves to Y -segments. 



SO 



V 



A 



V(/32-C^)(u2/3^C^-l) 
PC 



i 



(3.41) 



Obviously, it is A ^ C in general so (A, ^ (C, — i.e. the subsequent switching point 
is different from the symmetric image of the previous switching point with respect to 
the xi-axis. □ 

In the following proposition we use Lemma 13.71 to determine the form of the 
optimal trajectory. 

Proposition 3.8 (Form of the optimal trajectory). The optimal trajectory can 
have the one- switching form XY or the spiral form YX . . . YXY with an even number 
of switchings. 

Proof. We first show that when the optimal trajectory has more than one switch- 
ing, it cannot start with an X-segment. For just two switchings, consider the trajec- 
tory XYX depicted in Fig. 13.21 where a — 1 (starting point), (7, 0), 7 > 1 is the target 
point and (k,m),(CiC) the switching points. Since both of the switching points 
belong to the Y-segment passing through {(3,0), their coordinates satisfy (|3.18p . If 
we denote by s the common ratio 




then both k, C satisfy the equation 



(S + U2)xt - (M2/32 + -r^)xj -f 1 = 0, 



SO 




1 



S + U2 



< 1, 



since U2 > 1, s > 0. But also 




since > 1 and > 7^ > 1- Thus this trajectory cannot be optimal. 
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For more switchings, consider the case shown in Fig 13.11 where now a = I, and 
use s to denote the common ratio of the squares of the coordinates at the switching 
points. If T is the switching time between (C, ^) and (A, v), then from (j3.12p we obtain 

s cosh(2rY^) + 1 ^ 
ui cosh(2ry^) — 1 

But from l^^^ we find (a = 1) 

s + 1)(k2 _ 1) 



Ml UlK* 



< 1 4=^ (ui - 1)k^ > -1 



since ui > 1. Thus if the optimal trajectory has more than one switching, it needs to 
start with a K-segment. 

We next show that the optimal trajectory reaches the target point (7,0), 7 > 1 
with a y-segment. This is obviously the case for one switching, and also for two 
switchings since only the YXY trajectory is permitted (the XYX was excluded 
above) . For more than two switchings consider the situation shown in Fig. 13.21 It is 
fM^/K^ — ^ IC' = s and s > u\ since at least one FXy-segment is included in the 
trajectory. Point (C,0 belongs to the final X-segment ending to (7,0), so 

(s - + 72 = -wi7^ + -2 ■ 

The left hand side is positive, since s > ui, while the right had side is negative, since 
7 > 1, Ml > 1. Thus the optimal trajectory reaches the target point with a y-segment. 
□ 

Corollary 3.9. For \u\ < 1 the optimal solution has only one switching. 

Proof. For u — U2 — 1 the starting point (1,0) is an equilibrium point of system 
(|2.1ip . (j2.12|) . So the optimal trajectory cannot start with a y-segment. The only 
trajectory thus permitted is XY □ 

From Proposition 13.81 we see that the optimal trajectory can have aside from 



the expected one-switching form, shown in Fig. 3.3(a) the spiral form shown in Fig. 
3.3(b) An intuitive understanding of this latter form can be obtained by viewing 
system equations (|2.11[) . (|2.12p as describing the motion of a unit mass particle with 
position xi and velocity X2. In light of this interpretation we see that along a spiral 
trajectory the particle, instead of moving directly to the target, goes close to xi = 
where there is a strong repulsive potential {l/xf) to acquire speed and reach the target 
point faster. In the following theorem we calculate the transfer time for the candidate 
optimal trajectories. 

Theorem 3.10. Starting from (1,0), the necessary time to reach the target point 
(7,0),7 > 1 with one switching is 



Ml ly '-f^{ui + U2){ui + 1) J I y (mi + M2)(U27'' — 1) / 

(3.42) 

The necessary time to reach the target with n turns (2n .switchings) is 



Tn^Tj + nTx + {n- l)Ty + Tp, 



(3.43) 
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B(K,(i) 




(a) 




S2n('^2n'^'2n) 




(b) 



Fig. 3.3. (a) Trajectory with one switching (zero turns) (b) Trajectory with n turns. 



where 



T, = 



1 



2^ 
1 



: COS 



(s + U2) Vc? - 4U2 / ' 

-SCn+1 + U2^cl^+i - 4(s + U2) 



{s + U2)^l 



n+l 



AU2 



Tx 
Ty 



1 



1 



: cosh 



_1 / S + Ul 



27r — cos 



S — Ml 

I I S-U2 



S + U2 



Ci = M2 + 1, 



C„+l = M27 + 



(3.44) 
(3.45) 

(3.46) 
(3.47) 

(3.48) 
(3.49) 
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and s is the solution of the transcendental equation 



Cn+l + \/c2+l -4(S + W2) 



U2 



(3.50) 



in the interval ui < s < {u2 — 1)^/4. The constants ci and c„+i characterize the first 
and the last Y -segments, respectively, of the trajectory. The number of turns satisfies 
the following inequality 



n < 



To 



Tx{s+) 



(3.51) 



where s+ = {u2 



1)^/4 and [] denotes the integer part. 



Proof. In Fig. 3.3(a) we show a trajectory with one switching point B{k, fj,). The 

1 and /? = 7, 



coordinates of this point satisfy equations (|3.17|) and p.l8|) with a 
from which we find 

2 ^ M27'^ + i + 7^(m-l) 

72 (ui + U2) 



Using p.5|) with a = 1 and (13.61) with /3 = 7, we find that the necessary transfer 
time is given by (|3.42p . Next consider the case with n turns and 2n switching points 
{Kj,jjLj), Fig. 3.3(b)[ with constant ratio M^/^j = s. The first switching point satisfies 
the equations 



l4 



U2k\ 



UlK-l H 7j 



Cl, 



C, 



(3.52) 
(3.53) 



where ci is given by (|3.48p and c — —uia\ + 1/a^, while the second switching point 
satisfies 



lA 



M2 



U2n\ 



1 



C2, 



1 

2 ^' 
2 



(3.54) 
(3.55) 



where C2 = U2/3i + 1//??- The constants ci and C2 characterize the first and second Y- 
segments of the trajectory, while the constant c characterizes the X-segment joining 
them. Subtracting (|3.53l) from (I3.55|) and using Lemma lXTl which assures that ki 7^ K2 
(consecutive switching points are not symmetric with respect to xi-axis) we find that 



1 



s — Ui 



0. 



(3.56) 



But from p.52p . (|3.54p and the constant ratio relation we find 

2 



Cl + ^cf - 4(s + U2) ' 
C2 + \/c| - 4(s + U2) 

2(S + U2) ' 
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where, while solving the quadratic equations we used the — sign for the first and 
the + sign for the second switching point. The choice of sign for the first switching 
point will be justified below, while the choice of sign for consecutive switching points 
should be alternating to avoid picking the symmetric image of the previous point. 
Using these relations, (I3.56[) takes the form 

ci + \Jc{ - 4(s + U2) _ s - Ul 

C2 + \/4 - 4(S + U2) S + U2' 

By repeating the above procedure for all the consecutive pairs of switching points, we 
find 



, 6 — 1, Z, . . . , /t. 



U2 



Multiplying the above equations we obtain (13.501) . one transcendental equation for 
the ratio s. If we choose the + sign in the quadratic equation for k^, we obtain 
an equation similar to p.50p but with inverted left hand side. It is Cn+i > ci <^ 
(7^ - l)(u27^ - 1) > and ci, c„+i > 0, so 



c„+i + \l cl+^ - 4(s + U2) I s^u^ 

> 1 > ' 



Ci + ^Jcl-A{s + U2) \S + U2 

and the corresponding transcendental equation has no solution. Note that the left 
hand side of p.50p is a decreasing function of s while the right hand side is an 
increasing one, so if a solution exists, it is unique. The ratio is bounded below by the 
requirement s/ui > 1 and above by c^— 4(s + U2) > <^ s < s+ = (u2 — 1)^/4. This is 
also the maximum value of s on the first F-segment p.52p . Once we have calculated 
this ratio, we can find the time interval between consecutive switchings using (|3.46p for 
an X-segment and (|3.47p for a F-segment, relations obtained from Lemma on the 
inter-switching time. Observe that the times along all intermediate X- (respectively 
Y-) trajectories are equal. The initial time interval Tj (from the starting point up to 
the first switching) and the final time interval Tp (from the last switching up to the 
target point) can be easily calculated and are given in (|3.44l) and (I3.45p . respectively. 
The total duration r„ of the trajectory with n turns joining the points (1,0) and 
(7,0) is given by (I3.43p . Observe that T„(s) > nTx{s) > nTx{s+), where the last 
inequality follows from the fact that Tx is a decreasing function of s, see p.46p . A 
solution with n turns can be candidate for optimality only if the number of turns is 
bounded as in (j3.5ip . Otherwise we have r„(s) > Tq and the one-switching strategy 
is faster. □ 

We can find an approximate solution of (13.501) by setting s — ui (the lower limit) 
in the left hand side. We then obtain 

U1+U2VC 

Sn = — ;r7^: (3.57) 



1 - x/C 



where 



- ci + v^f^(Mi+M2) (3 53) 



Cn+l + ^Jcl+i -4(ui +U2) 
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Due to the monotonicity of the right and left hand sides of p.50|) . it is s„ > s„, 
where s„ is the exact solution. This approximation is good for s close to ui and thus 
for small n. As the xi coordinate 7 of the target point increases, the left hand side 
of ()3.50|) becomes less sensitive to variations in s, making the approximation more 
accurate. 

Using Theorem 13.101 we can find the times r„ for a specific target (7, 0) and 
compare them to obtain the minimum time. This is done in the next section for 
specific values of the control bounds. 



4. Examples. In Fig. 3.4(a) we plot the times To,Ti and T2 from Theorem 
13.101 corresponding to zero, one and two turns, for ui = 1,U2 = 8 and 76 [1. 10]. 
For 7 < 71 the strategy with zero turns (one switching) is optimal, while for 7 > 7i 
it is the strategy with one turn (up to the range of 7 plotted). The point (71, 0) can 
be reached with both strategies in equal time, that is, it belongs to the cut-locus [46] 
of these two control sequences from (1,0). Note that the strategies with one and 
two turns are feasible after some 7 > 1, where the transcendental equation p.50p 



has a solution. In Fig. 3.4(b) we plot the switching curves (black curves) as well as 



some characteristic optimal trajectories starting from (1,0). For 7 < 71 the optimal 
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trajectory starts with an X-segment that coincides with the switching curve (black 
curve) passing from (1,0). It switches at some point and then travels along a Y- 
segment (red curve) to meet the xi-axis. For 7 > 71 the optimal trajectory starts 
with a y-segment (red curve passing from (1,0)) and switches at some point in the 
tiny black area of this curve to an X-segment (blue curve). Then it meets at some 
point the second switching curve on the upper quadrant and changes to a F-segment 
(red curve) that hits the xi-axis at the target point. Note that the optimal trajectories 
between the two switchings (blue curves) are very close to the second switching curve 
on the upper quadrant and they are not shown entirely. 



In Fig. 3.5(a) we plot the times Tq, Ti, T2 and T3 from Theorem lS.lOl correspond- 
ing to zero, one, two and three turns, for ui — 1,U2 = 50 and 7 e [1,15]. Again, 
for small 7 the one-switching strategy is optimal and after some 7 = 71 the one-turn 
strategy becomes faster, but there is also a 7 = 72 beyond which the two-turn strat- 
egy is optimal (up to the range of 7 plotted) . The point (72 , 0) thus belongs to the 
cut-locus of the one- and two-turn control sequences from (1, 0) since it can be reached 
with one or two turns in equal time. In Fig. 3.5(b) we plot the switching curves (black 
curves) along with some characteristic optimal trajectories starting from (1,0). For 
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7 > 72 the optimal trajectory makes an additional turn. This is demonstrated by the 
three adjacent F-segments (red curves), which switch close to to the correspond- 
ing X-segments (blue curves), on a tiny switching curve which is hardly seen. In 
turn, these trajectories switch on the third switching curve on the upper quadrant to 
y-segments (red curves) that hit the xi-axis at the corresponding target points. 

5. Conclusion. In this article we formulated frictionless atom cooling in har- 
monic traps in an optimal control language and solved the corresponding time-optimal 
problem for a fixed initial condition (1, 0) and for varying terminal condition (7, 0), 7 > 
1. The optimal synthesis was obtained and an interesting switching structure was re- 
vealed. The results presented here can be immediately extended to the frictionless 
cooling of a two-dimensional Bose-Einstein condensate confined in a parabolic trap- 
ping potential [47] and even to the implementation of a quantum dynamical micro- 
scope, an engineered controlled expansion that allows to scale up an initial many-body 
state of an ultracold gas by a desired factor while preserving the quantum correlations 
of the initial state 05]. The above techniques are not restricted to atom cooling but 
are applicable to areas as diverse as adiabatic quantum computing and finite time 
thermodynamic processes. 

REFERENCES 

[1] H. MabuCHI and N. Khane,ia, Principles and applications of control in quantum systems, Int. 
J. Robust Nonlinear Control, 15 (2005), pp. 647-667. 

[2] A. C. DOHERTY, S. Habib, K. .Iacobs, H. Mabuchi, and S. M. Tan, Quantum feedback control 
and classical control theory, Phys. Rev. A, 62 (2000), 012105. 

[3] M. R. James, H. I. Nurdin, and I. R. Petersen, H°° control of linear quantum stochastic 
systems, IEEE Trans. Automat. Control, 53 (2008), pp. 1787-1803. 

[4] G. M. Huang, T. J. Tarn, and J. W. Clark, On the controllability of quantum-mechanical 
systems, J. Math. Phys., 24 (1983), pp. 2608-2618. 

[5] F. Albertini and D. D'Alessandro, Notions of controllability for quantum mechanical sys- 
tems, in Proceedings of the 40th IEEE Conference on Decision and Control, 2001, pp. 1589- 
1594. 

[6] C. Altafini, Controllability of quantum mechanical systems by root space decomposition of 
su(N), J. Math. Phys., 43 (2002), pp. 2051-2062. 

[7] A. Agrachev and T. Chambrion, Controllability and diameter of single-input quantum sys- 
tems, in Proceedings of the 1st IEEE Conference on Physics and Control, 3 (2003), pp. 720— 
725. 

[8] J.-S. Li and N. Khaneja, Ensemble control of Bloch equations, IEEE Trans. Automat. Control, 
54 (2009), pp. 528-536. 

[9] T. Chambrion, P. Mason, M. Sigalotti, and U. Boscain, Controllability of the discrete- 
spectrum Schroedinger equation driven by an external field, Annales de I'lnstitut Henri 
Poincare (C) Non Linear Analysis, 26 (2009), pp. 329-349. 

[10] A. M. Bloch, R. W. Brockett, and C. Rangan, Finite controllability of infinite- dimensional 
quantum systems, IEEE Trans. Automat. Control, 55 (2010), pp. 1797-1805. 

[11] K. BeauCHARD, J.-M. CorON, and p. RouOHON, Controllability issues for continuous- 
spectrum systems and ensemble controllability of Bloch equations, Commun. Math. Phys., 
296 (2010), pp. 525-557. 

[12] N. Khaneja, R. Brockett, and S. ,I. Glaser, Time optimal control in spin systems, Phys. 

Rev. A, 63 (2001), 032308. 
[13] D. D'Alessandro and M. Dahleh, Optimal control of two-level quantum systems, IEEE Trans. 

Automat. Control, 46 (2001), pp. 866-876. 
[14] U. Boscain, G. Charlot, J. -P. Gauthier, S. Guerin, and H.-R. Jauslin, Optimal control 

in laser-induced population transfer for two- and three-level quantum systems, J. Math. 

Phys., 43 (2002), pp. 2107-2132. 
[15] D. Stefanatos, N. Khaneja, and S. ,J. Glaser, Optimal control of coupled spins in the 

presence of longitudinal and transverse relaxation, Phys. Rev. A, 69 (2004), 022319. 



22 



D. STEFANATOS, H. SCHAETTLER, AND J.-S. LI 



[16] , S. J. Glaser, and N. Khaneja, Relaxation- optimized transfer of spin order in Ising 

spin chains, Phys. Rev. A, 72 (2005), 062320. 
[17] B. BONNARD AND D. SuGNY, Time-minimal control of dissipative two-level quantum systems: 

The integrable case, SIAM J. Control Optim., 48 (2009), pp. 1289-1308. 
[18] M. Lapert, Y. Zhang, M. Braun, S. J. Glaser, and D. Sugny, Singular extremals for the 

time-optimal control of dissipative spin 1/2 particles, Phys. Rev. Lett., 104 (2010), 083001. 
[19] D. Stefanatos and J.-S. Li, Constrained minimum- energy optimal control of the dissipative 

Bloch equations. Systems Control Lett., 59 (2010), pp. 601-607. 
[20] B. Bonnard, O. Cots, N. Shcherbakova, and D. Sugny, The energy minimization problem 

for two-level dissipative quantum systems, J. Math. Phys., 51 (2010), 092705. 
[21] A. Peirce, M. Dahleh, and H. Rabitz, Optimal control of quantum mechanical systems: 

Existence, numerical approximations, and applications, Phys. Rev. A, 37 (1988), pp. 4950— 

4964. 

[22] D. ,L Tannor, V. Kazakov, and V. Orlov, Control of photochemical branching: Novel 
procedures for finding optimal pulses and global upper bounds. Time Dependent Quantum 
Molecular Dynamics, Plenum, New York, 1992, pp. 347-360. 

[23] N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbruggen, and S. J. Glaser, Optimal 
control of coupled spin dynamics: Design of NMR pulse sequences by gradient ascent 
algorithms, J. Magn. Reson., 172 (2005), pp. 296-305. 

[24] ,J.-S. Ll .L Ruths, and D. Stefanatos, A pseudospectral method for optimal control of open 
quantum systems, J. Chem. Phys., 131 (2009), 164110. 

[25] T. Schulte-Herbruggen, S. J. Glaser, G. Dirr, and U. Helmke, Gradient flows for opti- 
mization in quantum information and quantum dynamics: Foundations and applications, 
Rev. Math. Phys., 22 (2010), pp. 597-667. 

[26] L I. Maximov, J. Salomon, G. Turinici, and N. C. Nielsen, A smoothing monotonic con- 
vergent optimal control algorithm for nuclear magnetic resonance pulse sequence design, 
J. Chem. Phys., 132 (2010), 084107. 

[27] J.-S. Li, J. Ruths, T.-Y. Yu, H. Arthanari, and G. Wagner, Optimal pulse design in 
quantum control: A unified computational method, Proc. Natl. Acad. Sci. U.S.A., 108 
(2011), pp. 1879-1884. 

[28] F. MOTZOi, J. M. Gambetta, S. T. Merkel, and F. K. Wilhelm, Coarse-grained optimal 
control methods for fast time-varying Hamiltonians, e-print arXiv: 1102. 0584 [quant-ph] . 

[29] T. Caneva, T. Calarco, and S. Montangero, Chopped random basis quantum optimization, 
e-print arXiv: 1103.0855 [quant-ph]. 

[30] C. E. WiEMAN, D. E. Pritchard, and D. J Wineland, Atom cooling, trapping, and quantum 
manipulation. Rev. Mod. Phys., 71 (1999), pp. S253-S262. 

[31] J. L Cirac and p. Zoller, New frontiers in quantum information with atoms and ions. 
Physics Today, 57 (2004), pp. 38-44. 

[32] H. R. Lewis and W. B. Riesenfeld, An exact quantum theory of the time- dependent harmonic 
oscillator and of a charged particle in a time- dependent electromagnetic field, J. Math. 
Phys., 10 (1969), pp. 1458-1473. 

[33] A. E. Leanhardt, T. A. Pasquini, M. Saba, A. Schirotzek, Y. Shin, D. Kielpinski, D. E. 

Pritchard, and W. Ketterle, Adiabatic and evaporative cooling of Bose-Einstein con- 
densates below 500 picokelvin. Science, 301 (2003), pp. 1513-1515. 

[34] S. BizE, P. Laurent, M. Abgrall, H. Marion, I. Maksimovic, L. Cacciapuoti, J. Grn- 

ERT, C. VlAN, F. PeREIRA DOS SANTOS, P. ROSENBUSCH, P. LeMONDE, G. SaNTARELLI, 

p. Wolf, A. Clairon, A. Luiten, M. Tobar, and C. Salomon, Cold atom clocks and 
applications, J. Phys. B: At. Mol. Opt. Phys., 38 (2005), pp. S449-S468. 
[35] D. Aharonov, W. van Dam, J. Kempe, Z. Landau, S. Lloyd, and O. Regev, Adiabatic 
quantum computation is equivalent to standard quantum computation, SIAM J. Comput., 
37 (2007), pp. 166-194. 

[36] P. Salamon, K. H. Hoffmann, Y. Rezek, and R. Kosloff, Maximum work in minimum time 
from a conservative quantum system, Phys. Chem. Chem. Phys., 11 (2009), pp. 1027-1032. 

[37] X. Chen, A. Ruschhaupt, S. Schmidt, A. del Campo, D. Guery-Odelin, and J. G. Muga, 
Fast optimal frictionless atom cooling in harmonic traps: Shortcut to adiabaticity, Phys. 
Rev. Lett., 104 (2010), 063002. 

[38] D. Stefanatos, ,J. Ruths, and J.-S. Li, Frictionless atom cooling in harmonic traps: A time- 
optimal approach, Phys. Rev. A, 82 (2010), 063422. 

[39] V. P. Ermakov, Second-order differential equations: Conditions of complete integrability, 
Appl. Anal. Discrete Math., 2 (2008), pp. 123-145, Translated from Russian by A.O. 
Harin, under redaction by P.G.L. Leach. 

[40] Y. Kagan, E. L. Surkov, and G. V. Shlyapnikov, Evolution of a Bose- condensed gas under 



TIME-OPTIMAL FRICTIONLESS ATOM COOLING 



23 



variations of the confining potential, Phys. Rev. A, 54 (1996), pp. R1753-R1756. 

[41] L. S. PONTRYAGIN, V. G. BOLTYANSKII, R. V. GAMKRELIDZE, AND E. F. MiSHCHENKO, The 
Mathematical Theory of Optimal Processes, Interscience Publishers, New York, 1962. 

[42] H. ,J. SuSSMANN, Time-optimal control in the plane, in Feedback Control of Linear and Nonlin- 
ear Systems, Lecture Notes in Control and Information Sciences, Vol. 39, Springer Verlag, 
Berlin, 1982, pp. 244-260. 

[43] H. J. SuSSMANN, The structure of time-optimal trajectories for single-input systems in the 
plane: The C°° nonsingular case, SIAM J. Control Optim., 25 (1987), pp. 433-465. 

[44] U. BOSCAIN AND B. PiCOOLl, Optimal Syntheses for Control Systems on 2-D Manifolds, 
Springer, SMAI, 2004. 

[45] V. JURDJEVIC, Geometric Control Theory, Cambridge University Press, Cambridge, 1997. 
[46] M. Berger, a Panoramic View of Riemannian Geometry, Springer- Verlag, Berlin, 2003. 
[47] ,1. G. MuGA, X. Chen, A. Ruschhaupt, and D. Guery-Odelin, Frictionless dynamics of 

Bose-Einstein condensates under fast trap variations, J. Phys. B: At. Mol. Opt. Phys., 42 

(2009), 241001. 

[48] A. del Campo, Frictionless quantum quenches in ultracold gases: A quantum dynamical mi- 
croscope, e-print farXiv: 1103.0714 [cond-mat. quant-gas]. 



